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ABSTRACT 



Aims. In this paper we use deep CFHT and SUBARU uBVRIz archival images of the Abell 383 cluster (z=0.187) to estimate its mass 
by weak lensing. 

Methods. To this end, we first use simulated images to check the accuracy provided by our KSB pipeline. Such simulations include 
both the STEP 1 and 2 simulations, and more realistic simulations of the distortion of galaxy shapes by a cluster with a Navarro- 
Frenk- White (NFW) profile. From such simulations we estimate the effect of noise on shear measurement and derive the correction 
terms. The /?-band image is used to derive the mass by fitting the observed tangential shear profile with a NFW mass profile. 
Photometric redshifts are computed from the uBVRIz catalogs. Different methods for the foreground/background galaxy selection are 
implemented, namely selection by magnitude, color and photometric redshifts, and results are compared. In particular, we developed 
a semi-automatic algorithm to select the foreground galaxies in the color-color diagram, based on observed colors. 
Results. Using color selection or photometric redshifts improves the correction of dilution from foreground galaxies: this leads to 
higher signals in the inner parts of the cluster We obtain a cluster mass Mvir = 7.5^^^ x lO'"* Mq. such value is ~ 20% higher than 
previous estimates, and is more consistent the mass expected from X-ray data. The ^-band luminosity function of the cluster is finally 
computed, giving a total luminosity L,„t = (2.14 ± 0.5) x 10'^ Lq and a mass to luminosity ratio M/L ~ 3OOM0/LO. 

Key words. Galaxies: clusters: individual: Abell 383 - Galaxies: fundamental parameters - Cosmology: dark matter 



1. Introduction 

Weak gravitational lensing is a unique technique that allows to 
probe the distribution of dark matter in the Universe. It mea- 
sures the very small distortions in the shapes of faint back- 
ground galaxies, due to foreground mass structures. This tech- 
nique requires a very accurate measurement of the shape pa- 
rameters as well as the removal of the systematic effects af- 
fecting them. In addition, galaxies to be used in the weak lens- 
ing analysis must be carefully selected so that they do not in- 
clude a significant fraction of unlensed sources, with redshift 
smaller than that of the lens. This would introduce a dilution 
and therefore an un derestimated signal, in particular towards the 
cluster center (see Broadhurst et al. 2005): such effect may be 
the reason of the under-prediction by weak lensing of the ob- 
served Einstein radi us, by a factor of ~ 2.5 (ISmith et al.ll200lt 
iBardeau et al.ll200"5l) . The optimal case would happen if photo- 
metric redshifts were available. Even if for weak lensing high 
accuracy in their estimate are not required for individual galax- 



* Based on: data collected at Subaru Telescope (University of Tokyo) 
and obtained from the SMOKA, which is operated by the Astronomy 
Data Center, National Astronomical Observatory of Japan; observa- 
tions obtained with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which 
is operated by the National Research Council (NRC) of Canada, the 
Institute National des Sciences de TUnivers of the Centre National de 
la Recherche Scientifique of France, and the University of Hawaii. This 
work is based in part on data products produced at TERAPIX and the 
Canadian Astronomy Data Centre as part of the Canada-France-Hawaii 
Telescope Legacy Survey, a collaborative project of NRC and CNRS. 



ies, on average we need at least 0-7/(1 + z) < 0.1: this implies 
having observations in several bands, spanning a good wave- 
length range. If few bands are available, an uncontaminated 
background sample can be obtained b y selecting only galaxie s 
redder than the cluster red sequence (iBroadhurst et al.ll2005h . 
However, such method often does not allow to get a number 
density of background sources high enough to allow an accu- 
rate weak lensing measure. In cluding galaxies bluer than the red 
sequence dOkabe et al.ll201C)l) requires a careful selection of the 
color offset, as bluer galaxies can be still contaminated by late- 
types memb ers of the cluster. Finally , if more than two bands 
are available. Medezinski et al.l(l2010h discussed how to identify 
cluster members and the foreground population as overdensities 
in the color-color space. 

In this paper we exploit deep uBVRIz images of the clus- 
ter Abell 383, taken with the MEG AC AM and SUPRIME cam- 
era mounted on the 3.6m CFHT and 8m SUBARU telescopes 
respectively, and publicly available. The mass of the cluster is 
derived by weak lensing, and values obtained by different selec- 
tion methods are compared. The properties of the cluster are re- 
viewed in Sect|2] Data reduction is discussed in Sect|3] In Sect.|4] 
we describe the algorithm used for the shape measurement, and 
some improvements for the removal of biases. The accuracy in 
the mass estimate is derived by comparison with simulations. In 
Sect. |5] we first summarize the different methods for the selec- 
tion of the background galaxies from which the lensing signal 
is measured. Such methods are applied to the case of Abell 383, 
and the masses derived in such way are then compared. Finally, 
in Sect. |6] we compare the mass derived in this paper with lit- 
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erature values, both by X-rays and weak lensing; a comparison 
is also done with the mass expected for the 7?-band luminosity, 
derived from the luminosity function of Abell 383. 

A standard cosmology was adopted in this paper: Qa - 0.7, 
= 0.3, Ho = 70 km s"' Mpc"', giving a scale of 2.92 
kpc/arcsec at the redshift of Abell 383. 

2. Abell 383 

Abell 383 is an apparently well relaxed cluster of gal axies of 
richness class 2 and of Baut z-Morgan type II-III ( Abell et alj 
I1989I) . located at z = 0.187 dFetisova et al.lll993h . It is domi- 
nated by the central cD galaxy, a blue-core emission-line Bright 
Cluster Gal a xy (B CG), that is aligned with the X-ray peak 
ISmith et all (12001 !). Abell 383 is one of the clusters of the 
XBACs sample (X-ray-Brightest Abell- type Clusters ), observed 
in the ROSAT All-Sky Survey (RASS: l\^gesl[T992h : its X-ray 
luminosity is 8.03 x lO"*'* erg sec ' in the 0.1 -2.4 keV band and 
its X-ray temperature is 7.5 keV (Ebelin g" et al.lll996l) . A small 
core radius, a steep surface brightness profile and an inverted de- 
projected temperature profile, show evidence of the presence of 
a cooling flow, as supporte d by the strong em ission Unes in the 
optical spectra of its BCG ( Rizza et al.]|1998l) . 

An extensive study of this cluster was carried out by 
ISmith et al] (1200 lb . in which lensing and X-ray properties were 
analyzed on deep optical HST images and ROSAT HRI data, 
respectively. A complex system of strong lensed features (a gi- 
ant arc, two radial arcs in the center and numerous arclets) were 
identified in its HST images, some of which are also visible in 
the deep SUBARU data used here. 

3. Data retrieval and reduction 

The cluster Abell 383 was observed with the SUPRIME cam- 
era mounted at the 8m SUBARU telescope: SUPRIME is 
a ten CCDs mosaic, with a 34 x 27 arcmin^ field of view 
(iMivazaki e t alJ l2002h . Data are publicly available in the BVRIz 
filters, with total exposure times of 7800s (R), ~ 6000s (B, V\ 
3600s (/) and 1500s (z); they were retrieved using the SMOKAQ 
Science Archive facility. Data were collected from seven differ- 
ent runs, amounting to ~ 55GB. Details about the observation 
nights and exposure times for each band are given in Table[T] The 
data reduction was done using the VST-Tube imaging pipeline, 
specific ally deve loped for the VLT Survey Telescope (VST, 
[Capacci oli et al.l l2005) but adaptable to other existing or future 
multi-CCD cameras (Grad o et al.ll20Ill) . 

The field of Abell 383 was also observed in the u* band with 
the MEGACAM camera attached on the Canada-France Hawaii 
Telescope (CFHT), with a total exposure time of 10541s. The 
preprocessed images were retrieved from the CADC archive. 

The basic reduction steps were performed for each frame, 
namely overscan correction, flat fielding, correction of the ge- 
ometric distortion due to the optics, and sky background sub- 
traction. It is worth mentioning that to improve the photometric 
accuracy, a sky superflat was used. 

Geometric distortions were first removed from each expo- 
sure, using the Scamp tool0 and taking the USNO-Bl as the as- 
trometric reference catalog. The internal accuracy provided by 
Scamp, as measured by positions of the same sources in differ- 
ent exposures, is ~ 0.05 arcsec. Different exposures were then 



' http://smoka.nao.ac.jp/ 

^ Stuff, Skymaker, SWarp and SExtractor are part of the 
Astromatic software developed by E. Berlin, see www.astromatic.net 
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Fig. 1. Density plots compaiing the photometric redshifts in the 
Abell 383 field available from the SDSS, and those here com- 
puted from the uBVRIz photometry. 
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Fig. 2. V - R vs R color plot: red and blue points are galaxies at 
Zphot = 0.187+ 0.1 and classified as early and late-type, respec- 
tively; dashed lines shows the +lcr levels. 



stacked together using SWarp. The coaddition was done in such 
a way that all the images had the same scale and size. 

The photometric calibration of the BVRI bands was per- 
formed using the standard Stetson fields, observed in the same 
nights of the data. In the case of the z band, we used a point- 
ing also covered by the Stripe 82 scans in the Sloan Digital Sky 
Survey (SDSS); the SDSS photometry of sources identified as 
point-like was used to derive the zero point in the SUBARU im- 
age. 

In the case of the MEGACAM-CFHT u* band images, re- 
duced images and photometric zero points are already available 
from the CADC public archive, hence only astrometric calibra- 
tion and stacking were required. 

Table |2] summarizes the photometric properties of the fi- 
nal coadded images (average FWHMs and limiting magnitudes 
for point-like sources) for each band. All magnitudes were 
converted to the AB system; magnitudes of sources classified 
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Table 1. Summary of obsei-vations done with the MEGACAM 
(m) and SUPRIME (BVRIz) cameras, used in this paper 



Date 


Band 


Exp. Time 


Total 


23 Dec 2003 


u' 


8381s 




21 Jan 2004 


u' 


2160s 


10541s 


09 Sept 2002 


B 


2400s 




08 Jan 2008 


B 


2400s 




09 Jan 2008 


B 


1200s 


6000s 


02 Oct 2002 


V 


4320s 




01 Oct 2005 


V 


1800s 


6120s 


12 Nov 2007 


R, 


2400s 




08 Jan 2008 


Rc 


5400s 


7800s 


02 Oct 2002 


Ic 


3600s 


3600s 


10 Sept 2002 


z' 


1500s 


1500s 



Table 2. Photometric properties of the coadded images; FWHM 
(arcsec) and Hmiting magnitudes (in the AB system) were com- 
puted for point-Hke sources. The signal to noise (SNR) is defined 
as SNR=FLUX.AUTO/FLUXERR^UTO. 



Band 


FWHM 


mag (SNR=5) 


mag (SNR=10) 


u* 


1.3 


26.1 


25.3 


B 


0.99 


27.0 


26.2 


V 


0.95 


26.6 


25.7 


Rc 


0.82 


27.1 


26.1 


Ic 


0.97 


25.0 


24.2 


z' 


0.74 


24.7 


23.9 



as galaxies we re corrected for Ga lactic extinction using the 
Schlegel maps (ISchlegel et al.lll998l) . 

The weak lensing analysis was done on the /?-band im- 
age. The masking of reflection haloes and diffraction spikes 
near bright stars are performed by ExAM, a code developed 
for this purpose. In short, ExAM takes the SExtractor cata- 
log as input, locates the stellar locus in the size-magnitude di- 
agram (see Sect. 14.1b . picks out stars with spike-like features 
from the isophotal shape analysis, and outputs mask region file 
that may be visualized in the DS9 software, and finally cre- 
ates a mask image in FITS format. The reflection haloes are 
masked by estimating the background contrast near the bright 
stars, whose positions are obtained from the USNO-B 1 . The ef- 
fective area available after removal of regions masked in such 
way was 801 arcmin^. Catalogs for the other bands were ex- 
tracted using SExtractor in dual-mode, where the /?-band im- 
age was used as the detection image. 

Photometric redshifts were c omputed from uBVRIz pho- 
tometry, using the ZEBRA code (iFeldmann et al.l l2006l) . This 
software allows to define 6 basic templates (elliptical, Sbc, 
Sbd, irregular and two starburst SEDs), and to compute log- 
interpolations between each pair of adjacent templates. We first 
applied to the BVz magnitudes the off'se t derived within the 
COSMOS survey by ICapak et all (l2007l) . fliat is h-0.19 (B), 
H-0.04 (V), -0.04 (z). We then convolved the stellar spectra from 
the Pickles' library (Pickles 1985) with the transmission curves 
used for each filter, and derived the offsets for the other filters, 
that is: 0.0 (m), 0.0 (R), +0.05 (I). Fig. [3] shows the comparison 




-0.5 0.5 1.5 -1.0 0.0 1.0 2.0 -0.5 0.0 0.5 1.0 

(B-V)ab (B-V)ab (V-R)ab 



Fig. 3. Observed (red dots) and model colors (black dots) for 
stars, after the offsets given in the text were applied. Model col- 
ors were derived convolving the Pickles' library of stellar spectra 
with the filter transmission curves. 
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Fig. 4. Distribution of the photometric redshifts computed from 
the uBVRIz data for R < 25 mag. 

of the model colors with those derived for the stars in our cata- 
logs, after the above offsets were applied. We also veiified that 
the offsets derived in this way are consistent with those obtained 
by running ZEBRA in the so-called photometry-check mode, 
that allows to compute magnitude offsets minimizing the aver- 
age residuals of observed versus template magnitudes. 

We then removed from the catalog those galaxies with a pho- 
tometiic redshift Zph > 3 and cr,/(l + z) > 0.1, where cr^ was 
derived from the 68%-IeveI errors computed in ZEBRA. The 
distribution of the so-obtained photometric redshifts is displayed 
in Fig. |4] The accuracy of the so obtained photometric redshifts 
was derived from the comparison with the SDSS photometric 
redshifts of galaxies in the SDSS (DR7), whose rms error is 
~ 0.025 for r' < 20 mag. We deiived (Fig. [TJ a systematic off- 
set Az/(1 + z) - 0.003 and an rms error cr Az/(1 + z) - 0.07. 
As a further check, we extracted from our catalog those galax- 
ies classified by ZEBRA as early-type, with R < 23 mag and 
kphot - 0.1 87 1 = 0.1. As expected, such galaxies define a red se- 
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quence (Fig|2]i: this was fitted as V - R - a + bR, with a - 0.5, 



4. Shape measurement 

Elhpticities of galaxies w ere estimated using the KSB approach 
(iLuppino & Kaise3ll997h : even if such algorithm does not allow 
to achieve accurate measurements of very low shear signals, y < 
10^^, it is nevertheless adeq uate in the case of w eak lensing by 
cluster s, as discussed e.g. bv lGill et al.l(l2009h and lRomano et al.l 
(IMS). 

In our KSB implementation, the SExtractor software was 
modified to compute all the relevant quantities, namely the raw 
ellipticity e, the smear polarizability P™, and the shear polariz- 
ability P^''. The centers of detected sources were measured using 
the windowed centroids in SExtractor 

The KSB approach assumes that the PSF can be described 
as the sum of an isotropic component (simulating the effect of 
seeing) and an anisotropic part. The correction of the observed 
ellipticity Cobs for the anisotropic part is computed as: 



^obs 



(1) 



where (starred terms indicate that they are derived from mea- 
surement of stars): 

P = <bs/^™*- (2) 

The intrinsic ellipticity e of a galaxy, and the reduced shear, g — 
7/(1 - A-), are then related by: 



(3) 



The term P^, introduced bv lLuppino & Kaiseil (Il997l) as the p re- 
seeing shear polarizability, describes the effect of seeing and is 
defined to be: 



psh* 



py _ psn _ psm^ — psh _ psm 



(4) 



The final output of the pipeline is the quantity eiso - 
e-dniso/P''', from which the average reduced shear, {g} = (eiso), 
provided that the average intrinsic ellipticity vanishes, (e) = 0. 

The ellipticity calculation is done by using a window func- 
tion in order to suppress the outer, noisy part of a galaxy: the 
function is usually chosen to be Gaussian with size 6. The size 
of the window function is commonly taken as the radius con- 
taining 50% of the total flux of the galaxy (as given by e.g. the 
FLUX_RADIUS parameter in SExtractor). In our case, we pro- 
ceed as follows. We define a set of bins with 6 varying between 2 
and 10 pixels (sources with smaller and larger sizes are rejected 
in our analysis), and a step of 0.5 pixel. For each bin we compute 
Cobs, P^^ and P™, and the elli pticity signal to noise ratio defined 
by Eq. 16 in lErben et al] jlOOll) : 



SNe(6') 



/ mw(\0\)d^e 



(5) 



The optimal size of the window function, ^max, is then defined as 
the value that maximizes SNe. Fig|5] shows the typical trend of 
SNe, which was normalized for display purposes, as a function 
of 0. It can be seen (Fig. |6]l that, on average, there is a con- 
stant off'set between 6^^^ and FLUX_RADIUS. Below SNe ~ 5, 
FLUX_RADIUS starts to decrease: this provides an estimate of 
limit on SNe, below which the shape measurement is not mean- 
ingful any more. 



6 [pixels] 



Fig. 5. SNe is displayed as a function of the window function 
size, 6, used to measure ellipticities. For display purposes, galax- 
ies were selected to have the same value of Omux, and SNe was 
normalized so that min(SNe)=0, max(SNe)=l. The vertical lines 
indicate the average (solid) and standard deviation (dashed) of 
FLUX_RADIUS for the same galaxies. 




Fig. 6. Running median of FLUXJlADIUS-^inax as a function 
of SNe. The vertical line shows the limit chosen for the selection 
of background galaxies. 



The terms p and q, derived from stars, must be evaluated 
at each galaxy posit ion: this is u sually done fitting them by a 
polynomial (see e.g. lRadovich et al. 2008), whose order must be 
chosen to fit the observed trend, without overfitting. The usage 
of the window function introduces a calibration factor, which 
is compensated by the P^ term. This implies that stellar terms 
must be co mputed and fi tted with the same value of 6 used for 
each galaxy (IHoekstra et al.l il998). An alternative approach, not 
based on a constant (and somehow arbitrary) order polynomial, 
is given e.g. by Generalized Additive Models: we found that 
the implementation in R (function GAM in the mgcv library) 
provides good results. Fig. |7] shows fitting and residuals of the 
anisotropic PSF component: from the comparison between the 
results obtained with polynomial and GAM fitting we see that 
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in the latter case we obtain lower residuals, in particular in the 
borders of the image. To quantify the improvement compared to 
the usage of the polynomial, we obtain (eaiiiso,i) = (1 ±5) x lO"'*, 
(eaniso2) - ("2 ± 9) X lO"'* with a polynomial of order 3, 
(eaniso!i> = (2 ± 4) X 10-4, (e,„i,o,2> = (1 + 7) X 10-4 with the 
GAM algorithm. The values of the fitted terms p and q at the po- 
sitions of the galaxies are predicted by GAM, that also provides 
an estimate of the standard errors of the predictions, Ap and Aq. 
From error propagation, the uncertainty on eiso was computed 
as: 

Ae? = (Ae,„i,o/P^)2 + (e3niso(P'')"'AP^)2, (6) 



where iAe,„i,,y = (Ae.^^y + (P™Ap)^ and {Apyy = (P' 
uncertainties on the measured values of f ™ and P^^ 
considered. 

For each galaxy, a weight is defined as: 

1 



Aqf- 
were not 



Aeo^ + Ael 



(7) 



where Aeo ~ 0.3 is the typical intrinsic rms of galaxy elliptici- 
ties. 



4. 1 . Star-galaxy classification 

Stars and galaxies were separated in the magnitude 
(MAG_AUTO) vs. size plot. Instead of using e.g. 
FLUX_RADIUS as the estimator of size, we used the quantity 
(5=MU_MAX-MAG_AUTO, where MU_MAX is the peak 
surface brightness above background. Saturated stars were 
found in the locus of sources with constant MU_MAX; in the 
6 vs. MAG_AUTO plot, stars are identified as sources in the 
vertical branch. Sources with 6 lower than stars were classified 
as spurious detections. In addition, we rejected those sources for 
which 6 is ~ 2cr larger than the median value. This is to exclude 
from the sample of stars used to compute the PSF correction 
terms, those sources for which the shape measurement may be 
wrong due to close blended sources, noise, etc. 

We further excluded those galaxies with w < 1 or SNe < 5, 
for which the ellipticity measurement is not meaningful. 

4.2. Error estimate in shear and mass measurement 

In the case of faint galaxies used for the weak lensing analy- 
sis, the ellipticity is underestimated due to noise. Such effect is 
not included in the P'' term, which c an be only com puted on 
stars with high signal to noise ratio. ISchrabback et al., (,2010i) 
proposed the following parametrization for such bias, as a func- 
tion of the signal to noise: 



= a * (SNe)*, 



(8) 



where and et are the ellipticities before and after the correc- 
tion respectively. Such parameters were derived using th e STEPl 
(iHevmans et al.il2006h and STEP2 (iMassev et alJ l200f) simula- 
tions, where both PSF and shear are constant for each simulated 
image. We obtain a - -0.1 and b - -0.45, which corresponds to 
a bias m changing from ~ 5% for SNe =5 to < 2% for SNe > 50. 
After such correction was applied, we computed again the aver- 
age shear from the STEPl and STEP2 simulations, and obtained 
a typical bias of ~ 3% for SNe = 5. 

We then estimated the accuracy on the mass that can be ob- 
tained from an image with the same noise and depth as in the 
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Fig. 7. PSF anisotropy correction derived with the GAM algo- 
rithm: the first three panels show the ellipticity pattern (mea- 
sured, fitted and residuals; X and Y are in pixels). The scale is 
displayed by the arrows in the upper right part of each panel 
(e - 0.05). In the next panel, black dots are the measured val- 
ues, green dots are after the correction; values rejected during 
fitting are marked in red. The last row shows for comparison the 
corrected ellipticities obtained using for the fit a polynomial of 
order 3. 



/?-band SUBARU image. To this end we dropped the assump- 
tion on constant shear, and produced more realistic simulations: 
the effect on galaxy shapes by weak lensing from a galaxy cluster 
was produced using the SHUFF code, that will be described in a 
separate paper (Huang et al., in preparation). To summarize, the 
code takes as input a catalog of galaxies produced by the Stuff 
tool; it computes the shear produced by a standard mass profile 
(e.g. Navarro-Frenk- White, NEW hereafter) and applies it to the 
ellipticities of the galaxies behind the cluster Such catalog is 
then used in the SkyMaker software, configured with the tele- 
scope parameters suitable for the SUBARU telescope and with 
the exposure time of the 7? -band image, producing a simulated 
image; the background rms of such image was set to be as close 
as possible to that of the real image. We considered for the lens 
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Table 3. Masses derived by simulations with NFW model fitting 
(M3d) and aperture densitometry (M2d)- Input masses used for 
each simulation are given in the first column (Min). 







M2D/M3D 




lO'^^Mo 




0.316 


0.31 ±0.13 


1.3 ±0.9 


1.00 


1.01 ±0.24 


1.4 ±0.4 


3.16 


3.12 ±0.36 


1.4 ±0.3 


10.0 


10.1 ±0.5 


1.3 ±0.1 



a range of masses at logMvir/M© - 13.5, 14, 14.5, 15.0, and a 
NFW mass profile with Cyir - 6. Each simulation was repeated 
50 times for each mass value, randomly changing the morphol- 
ogy, position and redshift of the galaxies. 

For each of these images, we run our lensing pipeline with 
the same configuration used for the real data. The density of 
background galaxies used for the lensing analysis was ~ 20 
gals arcmin~^. The fi t of t he mass was done as described in 
iRadovich et"al] (l2008h and iRomano et al] (|2010|) : the expres- 
si ons for the radial d epe ndence of tangential shear 77 derived 
bv lBartelman"nl (ll996') and' Wright & Brainerdl(l2000i) were used, 
and the NFW parameters (Mvir, Cvir) were derived using a max- 
imum likelihood approach. In addition, the 2D projected mass 
can be derived in a non-parametric way by aperture densitome- 
try, wher e the mass profile of the cluster is com puted by the ^ 
statistics (iFahlman et al ] ll994IClowe etalJll998h : 



^(0i) = ^(0<0i)-/c(02<0<0out) = 2 {yT)dln9 



{jT)dln0. 



(9) 



1 - (Ol/Oout^ 



f 



The mass is estimated as Afap(^i) — 7r^j^(^i)2crit, and ^out is 
chosen so that k(02. Bom) ~ 0. 

The average errors on mass estimate obtained in such way 
are displayed in Table |3] showing that masses can be estimated 
within an uncertainty of < 20% for M > 10'"* Mq. Such accuracy 
only includes the contribute due to shape measurement and mass 
fitting method, but it does not include the uncertainty due to the 
selection of the lensed galaxies. 

Finally, the masses derived by aperture densitometry are ~ 
1.3 higher than those o btained by mass fitting: this is in agree- 
ment with lOkabe et all qOlOl who find M2D/M30 = 1-34 for 
virial overdensity. 



5. Results 

One of the most critical source of systematic errors, which can 
lead to an underestimation of the true WL signal, is dilution of 
the distortion due to the contamination of the background galaxy 
catalog by unlensed foregroun d and cluster member galaxies 
(see e.g. lBroadhurst et aDl2005h . The dilution effect increases as 
the cluster-centric distance decreases because the number den- 
sity of cluster galaxies that contaminate the faint galaxy catalog 
is expected to roughly follow the underlying density profile of 
the cluster. Thus, correcting for the dilution effect is important 
to obtain unbiased, accurate constraints on the cluster parameters 
and mass profile. 

jBroadhurstetal.1 (l2005l) : lOkabe et al.l 



As discussed bv 

(l20Tcl) : IOguriiran 12010^ the selection of background galaxies 
to be used for the weak lensing analysis can be done taking those 



galaxies redder than the cluster red sequence. However, such 
selection produces a low number density (10 galaxies/arcmin^ 
in our case), and correspondingly high uncertainties in the de- 
rived parameters. In the following, we compare the results ob- 
tained by different methods. We first assumed that no informa- 
tion on the redshift is available, and that photometry from only 
one band is available (magnitude cut), or from more than two 
bands (color selection). Finally, we included in our analysis the 
photometric redshifts. The density of background galaxies is ~ 
25-30 galaxies/arcmin^, see Table|4] 

In order to derive the mass, we need to know the critical 
surface density: 



1 



47rG DiDis AjtG Dfi' 



(10) 



Dis, Ds, and D/ being the angular distances between lens and 
source, observer and source, and observer and lens respectively. 
This quantity should be computed for each lensed galaxy. As 
the reliability of photometric redshifts for the faint background 
galaxies is not well known, we prefer to adopt the single sheet 
approximation, where all background galaxies are assumed to 
lie at the same redshift, defined as /3(zs) = (J3(z)}- In the case 
of the selection based only on magnitude or colors, such value 
was derived from th e COSMOS catalog of photometric redshifts 
(ICapak et al.ll2007h . to which the same cuts used for the Abell 
383 catalog are applied. Later on, we computed /3{zs) from the 
photometric redshifts themselves, and compared such two val- 
ues. 

The mass was computed both by fitting a NFW profile 
(M3D - Mvir), and by aperture densitometry (M2d)- In the case 
of the NFW fit, in addition to a 2-parameters fit (virial mass Mvir 
and concentration Cvir), we also show (MNFW) the results ob- 
tained using the relation proposed by IBuUock et al.1 (1200 ih be- 
tween Mvir, Cvir and the cluster redshift Zc\- 



K IM, 



1 +Zc 



vir 

m7 



(11) 



withM* = 1.5 X IO'V/jMq, = 9, a = -0.13. 

The shear profiles obtained from the different methods dis- 
cussed below are displayed in Fig. [HJ also displayed are the av- 
erage values of tangential and radial components of shear, com- 
puted in bins selected to contain at least 200 galaxies, and cen- 
tered on the B CG: this is also wh ere the peak of the X-ray emis- 
sion is located dRizza et alj|1998l) . In order to check the possible 
error introduced by a wrong center, we considered a grid arou nd 
the position of the BCG, with a step of 2 arcsec: for each posi- 
tion in the grid, we took it as the center, performed the fit and 
derived the mass. Within 30 arcsec, we obtain that the rms is 
cr(Mvir) < 5%. The NFW parameters obtained by model fitting, 
and the reduced computed from the binned average tangential 
shear, are given in Table |4] 

The magnitude cut is the simplest approach as it only re- 
quires photometry from the same band in which the lensing mea- 
surement is done. Taking galaxies in the range 23 < 7? < 26 mag 
(a), produces a sample dominated by faint background galaxies, 
but the inner regions of the cluster may still present an unknown 
contamination by cluster galaxies. 

To improve the selection, we proceeded as follows. The locus 
of foreground galaxies was first found, allowing a better separa- 
tion of different galaxy populations (see'Medez inski et al.ll2010l 
and references therein), compared e.g. to methods based on the 
selection of only red galaxies. Here we considered two colors 
selections, namely B - z vs. V - z (b) and B - V vs. V - I (c). 
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Fig. 8. Shear profiles obtained with the different selection methods (see Table |4] where the parameters derived for each model are 
given). The fitting was done using the maximum likelihood approach. Binned points are shown for display purposes only. In panel 
a, the curves obtained by the two models (nfw/mnfw) overlap. 



with 21 < R < 26 mag. For intermediate redshifts (z ~ 0.2), 
foreground and background objects are well separated in such 
two colors diagrams. We explored the possibility to find the 
best selection criteria based only on the observed colors, with- 
out any information on the redshift of the galaxies. To this end, 
we developed a semi-automatic procedure, implemented in the 
R language dR Development Core Teaml2010l) . We first selected 
bright (R <21 mag) galaxies, which are expected to be mainly at 
Zph < 0.2. A kernel density estimate, obtained by the kde2d pack- 
age m R (IVenables & Riplevll2002h was then applied to these 
points, giving the plots displayed in Fig|9) the normalization was 
done so that the value of the maximum density in the binned data 
was equal to one. The boundary of the foreground galaxy region 
was then defined by the points within the same density level I 
(e.g., I = 0.2). Such region was converted to a polygon using the 
splanc^ package in R, that also allows to select for a given cata- 
log those sources whose colors lie inside or outside the polygon. 
A comparison with the model colors obtained in ZEBRA from 
the convolution of the spectral templates with filter transmission 
curves, shows that colors inside the area selected in such way 
are consistent with those expected for galaxies at redshift < 0.2. 
Galaxies classified as foreground in such way were therefore ex- 
cluded from the weak lensing analysis. 

We finally used photometric redshifts (d), both for the selec- 
tion of background galaxies, defined as those with 21 < R < 26 
mag, 0.3 < Zph < 3, and to compute the average value of /3: we 

^ http://www.maths.lancs.ac.uk/fowlings/Splancs 



obtain in such way f3(Zs) - 0.74, in good agreement with the 
value obtained from the COSMOS catalog with the same mag- 
nitude and redshift selection, j0(zj) - 0.73. 

The effect of the different selections on the residual pres- 
ence of cluster galaxies is displayed in Fig.IH] showing the den- 
sity of background galaxies computed in different annuli around 
the cluster: a clear increase of the density in the inner regions 
is visible in case a, which indicates that magnitude selection 
alone does not allow to completely remove the contamination by 
cluster galaxies. Such contamination is greatly reduced by color 
selection, and the optimal result is given by photometric red- 
shifts, as expected. As a further check, we also found for each 
method that the tangential shear signals of the rejected 'fore- 
ground/cluster' galaxies average out. In the following discus- 
sion, we take as reference the results from case d, which is very 
close to c in terms of uncertainties on fitted parameters, density 
of background galaxies and residuals in the radial component of 
the shear. 

It was pointed out in ' Hoekstral (l2003h and iHo ekstra et al.l 

that large-scale structures along the line of sight provide 
a source of uncertainty on cluster masses derived by weak lens- 
ing, which is usually ignored and increases as a larger radius 
(6'max) is used in the fitting. The uncertainty introduced by such 
component on the mass estimate can be ~ 10-20% for a cluster 
with M - IO^^Mq a t z ~ 0.2, g^ax - 10 arcmin as in our case 
(see Fig. 6 and 7 in iHoekstral (l2003b ). which is comparable to 
the uncertainties derived in the fitting. 
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Table 4. Best-fit NFW parameters: for each selection method, i n the upper row both Myir and Cvir were taken as free parameters, 
in the lower row the (Mvh, Cvir, Zci) relation ("Bulloc k et al.ll200"Tl) was used. The reduced was computed from the binned values 
displayed in Fig. [8] Also given are the values of M200: the larger uncertainties are due to the fact they were derived from the fitted 
values (Mvii , Cyir). The values of M2D before and after applying the factor 1 .34 (Sect J4.2t are displayed in the upper and lower rows, 
respectively. 



Method 


10'-* Mo 




^vir 

arcsec 


Mioo 
10'*Mo 


^7d.o.f. 


Mid 
10'*Mo 


Density 
gal/arcmin^ 


a 


7.783-3*21 


A nrilM 
^- '-'-1.49 




6.441-3 


1.13 


12.3 ±2.98 


30 




7.80if,2 


^- ' ^-0.12 


699\ 


6.46i-, 


1.13 


9.20 ± 2.22 




b 


'•"^--2.04 


5 772 33 


693**" 


6.42«f|, 


0.40 


8.14 ±3.15 


25 






4 62"- '3 


7331^,, 


7-44!f«6 


0.70 


6.07 ± 2.35 




c 


7 SI-'" 


5 40^ '3 


690!'6« 


6 30'' 

-3.08 


0.87 


10.5 ± 3.03 


28 




8.38if5, 




716«, 


6 942 "' 


0.61 


7.82 ± 2.26 




d 


7.54!f,, 


5-68!;'6o 


691.0513, 


6 ^S6 't5 
U.JJ_3Q2 


1.1 


10.2 ± 2.94 


25 




8.73if,„ 




725.77**47 


7 992.09 
;.^z,_l „ 


1.90 


7.59 ± 2.19 





Table 5. Luminosity function parameters, and uncertainties. <£* 
is in units of deg"^. 



Best-fit Standard error 

a TO6 007 

18.32 0.39 

0^/103 1.273 0.374 



„ . — < AO +0.55 

A weak l ensing analysis of Abell 383 was done by 
iBardeau et all (120071) using CFH12K data in the B, R, I filters. 
For the shape measurements they used a Bayesian method 
implemented into the IM2SHAPE software. To retrieve the 
weak lensing signal they selected the background galaxies 
as those within 21.6 < R < 24.9 mag and (R - /)> 0.7, 
obtaining a number density of ~ 10 gal arcmin"^. Their fit of 
the shear profile by a NFW profile gave as result a mass of 
M200 = 4.19 + 1.46 X /j^dlO'^Mo at rzoo = 1.32 + 0.17/i-,; Mpc 
and a concentration value of c = 2.62 + 0.69. 

Another weak lensing mass estimate of Abell 383 from 
CHF12K data was obtained by Hoekstra (2007) using two bands, 
B (7200 s) and R (4800 s). The background sample was se- 
lected by a magnitude cut 21 < R < 24.5, from which cluster 
red sequence galaxies were discarded. The remaining contam- 
ination was estimated from the stacking of several clusters by 
assuming that fraction of cluster galaxies fg^ was a function of 
radius oc r This function was used to correct the tang ential 
shear measurements. As discussed in lOkabe et al.l (12010 ^. this 
kind of calibration method does not allow to perform an unbi- 
ased cluster-by-cluster correction. Assuming a NFW profile, the 
fitted virial mass of Abell 383 was Mvir = 2.8if5 W^Mq. 



Also displayed in Table|4]are the projected masses computed 
by aperture densitometry from Eq. [TO] at a distance from the 
cluster center r - ryk, and 62 - 900", flout = 1000". A good 
agreement is found between these masses and the values com- 
puted from parametric fits, if we take into account the expected 
ratio M2D/M3D - 1.34 (see Sect.l4!2li. 

From the catalog based on photometric redshifts selec- 
tion (case d), we fi nally derived the S-map introduced by 
ISchirmer et~^ (|2004 . that is: S = M^p/crM.,^, where 

. ^-----'Q(l^' -^ol) (12) 



cr' 



i:,g>-6(i^--^oi) 
2(i:/w,)2 



The map was obtained by defining a grid of points along the 
image; the tangential components e,_, of the lensed galaxy ellip- 
ticities were computed taking as center each point in such grid. 
The weight w, was defined in Eq.|7] and 2 is a Gaussian function 
as in Radovic h et al. (2008): 

e(|,_,„|)._Lexp(-^^), (13) 

where 60 and 6s are the center and size of the aperture (0, ~ 
1.5 arcmin). The S-map is displayed in Fig.fTO] showing a quite 
circular mass distribution centered on the BCG. 



6. Discussion 

Several mass measurements of this cluster are available in liter- 
ature, based on different data and/or methods. ' Schmidt & AllenI 
(I2007h used Chandra data and modeled the dark mat- 
ter halo by a generalized NFW profile, obtaining a mass 
value Mvii- = 9.16;^} ?? x 10'"* and a concentration value 
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v-z 




o.e 

B-V 



Fig. 9. Color selection of foreground galaxies (Zphot < 0.2). The 
contour in red is the density level chosen for the selection; the 
points display the model colors computed by ZEBRA for this 
redshift range. 



Abell 383 belongs to the clusters sample selected for the 
Local Cluster Substructure Survey (LoCuSS) project (P.I. G. 
Smith). Within this project, a weak lensing analysi s of this clus- 
ter has been recently performed by Ok abe et al.l (j2010) using 
SUBARU data in two filters, /' (36 minutes) and V (30 min- 
utes). In addition to a magnitude cut 22 < i < 26 mag, they used 
the color information to select galaxies redder and bluer than the 
cluster red sequence. Looking at the trend of the lensing signal 
as function of the color offset, they selected the sample where 
the dilution were minimized, obtaining a background sample of 
~ 34 gal arcmin"^ for the computation of tangential shear pro- 
file of the cluster The fit of this profile by a NFW model did 
not give an acceptable fit: the virial mass computed assuming 
a NFW profile was Mvh- = 3.62+' '^ x h-^lO^^M^ with a high 
concentration parameter Cvii- - S.S?;^^ ^^. The same authors also 
derived the projected mass, obtaining M2D - 8.69 x /i"' lO'^'M© 
at the virial radius. 

Scaling such mass values to h - 0.7, we derive Mvir ~ 
(4 - 5) X 10'"^ Mq. This value is still consistent within uncer- 
tainties with the value derived in the present analysis (Mvii - 
(l.S^Jg X 10'"* Mq, Cvir = 5.7^j g, case d). In our case, we obtain 




41.99 42.00 42.01 42.02 42.03 42.04 



Fig. 10. Weak lensing S-map showing the mass distribution de- 
rived by weak lensing; overlaid is the central region of the Abell 
383 field. 



-1- 

•i- 



-5- 

-I- 



4 5 

r (arcmin) 



Fig. 11. Density (gal/arcmin^) of background galaxies used for 
the lensing analysis, as a function of the distance from the clus- 
ter, for the different selection methods here considered. The case 
with no selection is also displayed for comparison. 



however a better agreement with both the values of Myir and Cyir 
given by the X-ray data, and a better consistency between para- 
metric and non-parametric mass estimates, once the projection 
factor is taken into account. This may be due either to how the 
selection of background/foreground galaxies was done, or to a 
higher accuracy in the shape measurement as a combination of 
the deeper image and calibration of the bias due to SNR (Sect.|4]i. 
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Fig. 12. /?-band LF of the galaxies in Abell 383. Data points are 
derived by binning the data in magnitude bins of 0.5 mag and 
error bars by Poissonian errors; the curve traces the best-fitting 
LF and the shaded area marks the model uncertainty, obtained 
by a bootstrap technique. 



Finally, we computed the luminosity function (LF hereafter) 
in the /?-band and derived the total luminosity by integrating 
the fitted Sche chter (Schechter 1976) function as described in 
iRadovich et al.l feOOSh . In order to obtain the cluster LF, that is 
the number of galaxies per unit luminosity and volume belong- 
ing to the cluster, we need to remove from our catalog all the 
background and foreground galaxies. Usually, this is done by 
statistically subtracting the galaxy counts in a control field from 
galaxy counts in the cluster direction. Here we take advantage 
of the selection on the color-color diagrams described in Sect. |5] 
and extract a catalog that includes cluster, foreground and resid- 
ual background galaxies. The last two components were further 
removed by the statistical subtraction, where we defined as clus- 
ter area the circular region around the cluster center of radius 
r=10.2 arcmin (1.3Mpc); the control field was instead defined as 
the area outside the circle of radius r=15.3 arcmin. For the best- 
fitting procedure to the Schechter function, we adopted conven- 
tional routines of minimization on the binned distributions. Best- 
fitting parameters are listed in Table |5] The /?-band total lumi- 
nosity, calculated as the Schechter integral, is Ltot - (2.14 + 0.5) 
X 10'^ Lq. The errors were estimated by the propagation of the 
68%-confidence-errors of each parameter 

Fo r comparison, we then used the relation in iPopesso et al.l 
(l2007h between M200 and the optical luminosity. Lop, to see 
whether the mass obtained in this paper is consistent with the 
value expected for that luminosity. According to this relation, 
the mass expected for such luminosity is M20Q = (4.73 + 1.3)x 
lO''* Mq, in good agreement with the value derived by our weak 
lensing analysis, M200 ~ 6.3x lO'"* Mq (case d), corresponding 
to M/L ~ 3OOM0/LO. 



7. Conclusions 

We have computed the cluster Abell 383 mass by weak lens- 
ing, using a deep /?-band image taken with the Suprime cam- 
era on the Subaru telescope. Catalogs extracted from combined 



CFHT-hSUBARU uBVRIz images were used to derive photo- 
metric redshifts, and improve the weak lensing analysis. The 
data were reduced using a pipeline developed in-house, that was 
specifically designed for wide-field imaging data. The elliptic- 
ities, from which the shear signal was derived, were measured 
using a pipeline based on the KSB approach. We discussed some 
aspects that may improve the results, namely the size of the win- 
dow used to suppress the noise from the outer part of the galax- 
ies, the selection of a limit on SNR below which the measure- 
ment ellipticity is not accurate enough, and a weighting scheme 
where uncertainties on spatial fitting of the PSF correction terms 
were taken into account. 

The accuracy on the mass estimate by weak lensing available 
with our KSB pipeline was first derived on simulated images 
which were built in such a way to mimic as closely as possi- 
ble the background noise and the depth of the real image. From 
these simulations we conclude that the mass can be measured 
with an uncertainty ~ 5-10% for log M/Mq > 14.5. Such ac- 
curacy takes into account the measurement errors of the ellip- 
ticity, but not the errors due to e.g. the foreground/background 
galaxy separation, that may introduce an underestimate of the 
mass. The impact of such selection was evaluated by comparing 
three methods for the foreground/background galaxy separation, 
namely: magnitude cut in one band, color selection and usage of 
photometric redshifts. All methods gave consistent estimates of 
the total virial mass, but from the shear profile it can be seen that 
a dilution of the signal in the inner regions is still present, in the 
case of a simple magnitude cut. Color selection and photometric 
redshifts provide better results, even if the accuracy of the pho- 
tometric redshifts is not high due to the few available bands. The 
virial mass of Abell 383 here obtained by NFW model fitting is 
in agreement with the value obtained from the non-parametric 
mass estimate, that is Myir ~ 7 x 10'"* Mq. Other previous weak 
lensing analyses give Myir ~ (4-5) x 10'"^ Mq: the value found in 
this paper seems more in agreement with the value found by X- 
ray data, and we also have a better agreement betw een paramet- 
ric and non-parametric estimates, compared to e.g. lOkabe et aP 

Finally, we estimated the /?-band LF of Abell 383, and de- 
rived the total /?-band luminosity of the cluster: starting from this 
value and usin g the relation between mass and luminosity found 
for clusters bv lPopesso et al.l (l2007l) . we conclude that the mass 
derived by weak lensing is consistent with the value expected for 
this luminosity. 
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